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A GRAPHICS  PROGRAM  FOR  UPDATING 


THE  CONFIDENCE  REGION  OF  A TARGET 

INTRODUCTION 

The  capability  to  update  the  confidence  region  associated  with  a 
target  has  many  tactical  applications.  Useful  analytic  solutions  to 
this  problem  are  very  difficult  to  obtain  (see  for  example  [l]  and  [2]). 
The  graphics  program  presented  here  performs  the  dynamic  update  under 
fairly  general  conditions.  The  algorithms  employed  are  mathematically 
rigorous;  the  main  approximations  made  are  associated  with  replacing 
the  continuous  boundary  curve  of  a region  with  a discrete  set  of 
equally  spaced  points  on  the  curve. 

The  operations  of  the  Tenth  Fleet  during  World  War  II  has  shown 
that  information  collected  and  processed  by  a wide  area  ocean  surveil- 
lance system  can  have  significant  tactical  applications.  Human  judg- 
ment plays  a key  role  in  the  processing  of  tactical  ocean  surveillance 
data.  A graphics  terminal  is  a natural  interface  for  a man-machine 
tactical  information  processing  system.  Computers  have  been  used  with- 
in wide  area  ocean  surveillance  systems  for  the  direct  recall  of  infor- 
mation or  for  elementary  data  correlation,  but  their  capabilities  have 
not  been  fully  exploited  in  this  application.  A practical  tactical 
information  processing  system  should  satisfy  the  following  criteria: 

Note:  Manuscript  submitted  November  10,  1976. 


(a)  It  should  be  able  to  process  multi-source  information. 

(b)  It  should  be  an  interactive  system  which  can  be  used  and 
understood  by  non  mathematically  oriented  Navy  analysts. 

(c)  It  should  enable  the  analyst  to  interrelate  his  under- 
standing of  the  tactical  situation  with  the  performance 
characteristics  of  the  sensor  systems. 

The  confidence  region  update  procedure  which  is  described  in  this  report 
satisfies  the  above  requirements. 

DESCRIPTION  OF  UPDAT  PROGRAM 

The  program  UPDAT,  which  updates  the  confidence  region  associated 
with  a target  is  written  in  FORTRAN/PLOT  10  for  the  PDP-10  computer.  It 
updates  a confidence  region  associated  with  one  or  more  sensor  detects 
subject  to  specified  velocity  constraints.  The  three  basic  inputs 
required  are: 

(a)  An  initial  confidence  region  P in  position  space  which 
is  assumed  to  contain  the  target. 

(b)  A constraint  set  [ / in  velocity  space. 

(c)  An  update  time  T. 

The  UPDAT  program  computes  the  updated  confidence  region  P^  in 
position  space  after  an  elapsed  time  T subject  to  the  condition  that 
the  velocity  of  the  target  lies  in  the  constraint  set  1/  . The  bound- 
ary of  the  updated  region  is  displayed  on  a graphics  system. 

The  constraint  set  (/  is  always  defined  by  specifying  the  lower 
and  upper  bounds  on  speed  (VI  and  V2)  and  the  lower  and  upper  bounds 
on  heading  (A1  and  A2).  In  general  VI  < V2  and  A1  < A2.  Degenerate 
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cases  in  which  VI  = V2  and  A1  = A2  ( 1/  contains  a single  velocity 
vector,  i.e.,  there  is  no  uncertainty  in  velocity)  or  VI  = V2  and 
Al  < A2  (1/  contains  all  velocity  vectors  lying  on  a segment  of  a 
circle  centered  at  the  origin  in  velocity  space)  or  Al  = A2,  VI  < V2 
( 1/  contains  all  velocity  vectors  on  a segment  of  a ray  emerging  from 
the  origin  in  velocity  space)  are  all  permissible. 

The  analyst  operating  the  UPDAT  program  has  the  following  three 
options  for  the  selection  of  an  initial  confidence  region  P in  position 
space: 

(a)  Elliptical  confidence  region.  The  analyst  inputs  five 
parameters  which  define  an  arbitrary  ellipse. 

(b)  An  arbitrary  convex  confidence  region.  The  analyst  inputs 
the  boundary  points  of  this  set. 

(c)  A wedge  shaped  confidence  region.  The  analyst  inputs  four 
parameters  which  describes  an  arbitrary  wedge. 

OPERATIONAL  INPUTS  TO  UPDAT 

The  operational  inputs  to  UPDAT  are  either  of  the  sensor  perform- 
ance type  or  the  tactical  type.  The  performance  type  of  input  arises 
from  the  known  physical  characteristics  of  the  ocean  surveillance  (OS) 
sensors.  An  OS  sensor  measurement  of  a kinematic  variable  (e.g.,  such 
as  speed,  heading,  position,  line  of  bearing)  yields  an  expected  value 
of  this  variable  and  a confidence  region  (at  a selected  level)  about  the 
expected  value.  The  tactical  type  of  input  arises  from  the  analyst's 
understanding  of  the  tactical  context  of  the  situation. 

The  sensor  performance  inputs  arise  from  the  following  broad  classes 
of  sensors : 
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(a)  Position  measuring  sensors. 


(t)  Line  of  tearing  (LOB)  measuring  sensors. 

(c)  Velocity  measuring  sensors. 

Elements  within  the  first  class  include  HF/DF  and  SOSUS.  The  out- 
put of  this  class  of  sensors  is  typically  an  elliptical  confidence 
region.  TJPDAT  requires  the  following  five  input  parameters  to  describe 
an  elliptical  confidence  region 

(XO,  YO,  Rl,  R2,  SO) 

where 

XO  = longitude  of  center  of  ellipse, 

YO  = latitude  of  center  of  ellipse, 

Rl  = length  of  semi -major  axis, 

R2  = length  of  semi-minor  axis  (where  R^  > Rg  > 0),  and 
SO  = orientation  of  semi-major  axis  relative  to  north. 

The  second  class  of  inputs  includes  data  from  many  of  the  passive 
sensors  which  measure  LOB  information.  Elements  of  this  class  arise 
from  isolated  detects  at  individual  HF/DF  or  SOSUS  sites  as  well  as 
from  passive  SW  sensors.  There  are  four  input  parameters  which  can  be 
associated  with  a LOB  message: 

(XB,  YB,  31,  B2) 

These  parameters  determine  a wedge  shaped  confidence  region  as  follows: 

XB  = longitude  of  sensor  location 
YB  = latitude  of  sensor  location 

B1  = line  of  bearing  of  first  great  circle  of  wedge 

B2  = line  of  bearing  of  second  great  circle  of  wedge,  with  B2  > Bl. 

. 

The  coordinates  of  the  vertex  of  the  wedge  are  (.XB,  YB).  Let  30  denote 
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the  observed  bearing  and  let  A3  denote  the  maximum  bearing  error. 


Then 


B1  = BO  - A3 
32  = 3C  + aB 

The  determination  of  velocity  from  the  observed  Doppler  shift  of 
acoustic  frequency  is  one  source  of  velocity  information.  Velocity 
information  can  also  arise  from  visual  or  radar  detects.  Intelligence 
reports  (e.g.,  intercepted  communications)  provide  another  source  of 
velocity  data.  In  UPDAT  the  velocity  data  are  assumed  to  be  described 
by  the  four  quantities: 

(Al,  A2,  VI,  V2) 

where 

Al  = lower  bound  of  heading 
A2  = upper  bound  of  heading 
VI  = lower  bound  of  speed 
V2  = upper  bound  of  speed. 

With 

A2  > Al 
V2  > VI. 

Tactical  considerations  will  often  determine  the  constraint  set 
V = (Al,  A2,  VI,  V2) . Selected  bounds  (VI,  V2)  of  the  target's  speed 
could  originate  from  a combination  of  platform  performance  data  and 
tactical  context.  As  an  example,  the  analyst  should  have  available 
information  on  the  top  speed  WAX  and  the  cavitation  speed  VC  of 
different  submarine  classes.  Under  certain  tactical  situations  V2 
could  equal  VMAX  while  under  other  situations  V2  will  be  equal  to 
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VC.  Selected  bounds  (Al,  A2)  of  the  target's  heading  could  arise  from  a 
knowledge  of  the  target's  likely  objectives  or  its  observed 
behavior.  The  analyst  can  select  Al  = 0°  and  A2  = 36O0  if  there 
is  no  information  to  narrow  down  the  heading  estimates. 

•TACTICAL  APPLICATIONS 

Two  applications  of  the  UPDAT  system  at  a central  OS  site 
are  described  here.  In  the  first  application,  the  analyst  will 
use  the  UPDAT  output  to  assist  him  in  making  track  assignment  decisions. 
In  the  second  application,  the  UPDAT  system  is  used  to  develop  predic- 
tions of  the  locations  of  enemy  combatants, which  are  then  sent  to 
search  and  attack  units  (SAUs).  The  track  assignment  application  of  the 
UPDAT  system  would  develop  with  the  following  sequence: 

(a)  At  time  T an  OS  sensor  has  "detected"  a target.  This 

o 

contact  is  at  time  T the  last  reoortea  observation 

o 

of  the  target.  (The  target's  track  could  at  time  Tq 
consist  of  this  single  report.) 

(b)  At  time  T + T a second  detect  was  made.  It  is 

o 

assumed  that  the  signature  information  in  the  second 
detect  is  not  by  itself  sufficient  to  assign  it  to  the 
first  target. 

(c)  Using  the  procedure  outlined  in  the  proceeding  section 

a confidence  region  P in  position  space  and  a confidence 

region  \j  in  velocity  space  associated  with  the  target 

at  time  Tc  are  developed.  Let  ?'  denote  the 

elliptical  confidence  region  associated  with  the  new 

sensor  detect  at  time  T + T.  Let  P denote  the  uodate 

o u 
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to  time  Tq  + I of  the  set  p subject  to  the  velocity 
constraint  set  1/  . If  p/  doesn't  intersect  P 
then  the  assignment  of  the  second  report  to  the  existing 
track  would  not  be  consistent  with  the  analyst's  know- 
ledge of  the  target's  behavior.  If  P/  does  intersect 

p , the  assignment  is  consistent  with  his  knowledge 

u 

but  the  analyst  is  not  forced  to  make  this  decision. 
These  confidence  regions  are  sketched  in  the  figure  shown  below. 


POSITION  SPACE  VELOCITY  SPACE 

Fig.  1 — Elements  of  track  assignment  application 


The  tactical  prediction  application  would  proceed  as  follows:  At 

tine  T a target  detect  was  obtained.  As  discussed  in  the  preceeding 
paragraph  the  confidence  regions  P and  1/  are  developed.  SAUs  are  to 
attempt  to  acquire  the  target  at  time  Tq  + T.  The  UPDAT  system  can  be 
used  to  combine  the  confidence  regions  P and  \J  with  the  update  time 
T to  compute  the  updated  confidence  region  which  is  the  region 

where  the  search  effort  should  be  concentrated.  Given  the  shape  and 
size  of  and  the  sensor  and  platform  characteristics  of  the  SAUs, 

composite  search  strategies  can  be  evaluated  according  to  their 
associated  total  probability  of  detection. 

UPDAT  FLOWCHART 

The  flow  chart  shown  in  Figure  2 describes  the  flow  of  operations 
in  UPDAT.  The  input  set  (Al,  A2,  VI,  V2,  T)  determine  the  heading  and 
speed  bounds  of  the  target  and  the  update  time.  The  second  input  set 
(i)  determines  which  of  the  three  options  for  the  confidence  region  is 
selected. 

For  the  ellipse  case  the  next  input  set  (XO,  YO,  Rl,  R2,  SO) 
contains  the  coordinates  of  the  center  of  the  ellipse,  the  semi -major 
and  semi-minor  axis  and  the  orientation  of  the  ellipse.  The  subroutine 
30E  generates  30  points  lying  on  the  boundary  of  the  selected  ellipse. 
The  call  to  CONDRAW  plots  this  boundary  on  the  graphics  terminal.  The 
subroutine  TEST  determines  which  subroutine  (algorithm)  is  used  to 
perform  the  update  operation. 

If  the  uncertainty  in  position  space  is  zero  (R^  = Rp  = 0)  or  is 
small  compared  to  T times  the  uncertainty  in  velocity  space,  then 


UPDATE 
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subroutine  UPBA3  is  called.  This  subroutine  calls  subroutine  BOS  which 
in  turn  generates  the  boundary  points  of  the  constraint  set  in  velocity 
space.  UP BA 3 now  combines  this  set  of  velocity  points  with  the  single 
position  point  and  the  update  time  to  generate  the  updated  confidence 
region.  The  call  to  CONDRAW  plots  the  updated  confidence  region. 

If  the  uncertainty  in  velocity  space  is  zero  (A1  = A2  and  VI  = V2) 
or  is  small  compared  to  the  uncertainty  in  position  space  divided  by  T 
then  subroutine  UPDA4  is  called.  Subroutine  UPHA4  performs  a translation 
of  the  initial  region  in  position  space  as  determined  by  the  single 
velocity  vector  and  the  update  time . The  call  to  CONDRAW  plots  the  up- 
dated confidence  region. 

If  the  lower  bound  of  speed,  VI,  is  zero,  then  subroutine  UPBAl  is 
called.  For  VI  = 0,  the  constraint  set  in  velocity  space  will  be 
convex.  It  is  shown  in  Appendix  3 that  a convex  set  in  position  space 
when  combined  with  a convex  set  in  velocity  space  will  yield  an  updated 
sat  in  position  space  which  is  also  convex.  The  subroutine  BOP  first 
computes  an  updated  set  of  position  points  by  combining  all  points  on 
the  position  boundary  (determined  by  30E)  with  all  points  on  the 
velocity  boundary  (determined  by  BOS).  A subset  of  this  set  contains 
the  boundary  points  of  the  updated  convex  set.  This  set  of  boundary 
points  is  determined  using  the  algorithm  described  in  Appendix  C. 

For  the  general  case  (AR  ^ 0,  AV  £ 0,  VI  l 0)  subroutine  PART  is 
called.  Subroutine  PART  decomposes  V into  a union  of  almost  convex 
sets  as  determined  by  a preselected  convexity  criterion.  For  each  one 


of  these  nearly  convex  sets  UPDA1  is  called. 

The  case  where  1 = 2 is  very  similar  to  the  1=1  case.  Here 
the  boundary  points  of  the  initial  convex  set  are  entered  by  the  analyst. 
Subroutine  ARHIP  plays  nearly  the  same  role  as  subroutine  TEST.  The 
first  call  to  CONDRAW  plots  the  initial  confidence  region. 

For  the  wedge  case  the  next  input  set  (XB,  YB,  Bl,  32)  contains 
the  coordinates  of  the  vertex  of  the  wedge  and  the  two  angles  which 
define  the  great  circle  boundaries  of  the  wedge.  Except  for  near  the 
vertex  point  the  updated  confidence  region  will  be  bounded  by  two  great 
circles.  These  great  circles  are  determined  by  associating  with  the 
lines  of  bearings  Bl  and  32  two  critical  velocity  vectors  lying  in  p. 
The  two  critical  velocities  are  determined  exactly  and  hence  except  for 
near  the  vertex  point  WEDUP  computes  an  exact  update  of  the  initial 
wedge.  The  first  call  to  CONDRAW  plots  the  initial  confidence  wedge 
and  the  second  call  to  CONDRAW  plots  the  updated  confidence  region. 
EXECUTION  OF  UPDAT 

This  section  describes  how  a Navy  analyst  would  use  UPDAT.  The 
UPDAT  program  operates  in  an  interactive  mode.  Data  are  entered  on  line 
on  a graphics  terminal  with  the  program  providing  cues  for  the  required 
input  at  each  stage. 

The  system  requests  the  first  input  parameter  by  printing  out  the 
phrase 

HEIGHT: 

The  analyst  then  inputs  the  value  of  the  height  above  sea  level  in 
nautical  miles  of  the  point  of  observation. 

The  system  requests  the  second  set  of  input  parameters  by 
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printing  out  the  phrase: 

LAT,  LCN: 

The  analyst  then  inputs  'the  values  of  the  latitude  and  longitude  in 
(decimal)  degrees  of  the  point  of  observation. 

The  system  requests  the  third  set  of  input  parameters  by  printing 
out  the  phrase : 

BIFUT  Al,  A2,  VI,  V2,  T: 

The  analyst  then  inputs  the  five  quantities 

Al  = lower  bound  of  heading  of  target  (in  degrees) 

A2  = upper  bound  of  heading  of  target  (in  degrees 

VI  = lower  bound  of  speed  of  target  (in  knots) 

V2  = upper  bound  of  speed  of  target  (in  knots) 

T = update  time  (in  hours) 

The  heading  bounds  Al  and  A2  are  specified  clockwise  relative  to 
due  north.  If  the  spread  of  ’uncertainty  (Al,  A2)  in  target  heading 
does  not  span  due  north  then  Al  < A2.  If  the  spread  of  headings  does 
contain  due  north  then  Al  > A2. 

» The  system  now  requests  the  fourth  set  of  input  parameters  by 

printing  out  the  phrase : 

INPUT  I:  (1=1  FOR  ELLIPSE:  1=2  FOR  ANALYST  SELECTED  SET: 

I = 3 FOR  WEDGE ) 

The  analyst  then  inputs  the  value  of  I equal  to  1,  2,  or  3 depending 
on  which  option  he  selects  for  the  initial  confidence  region  of  the 
target ; 

If  the  analyst  inputs  the  value  of  1=1,  then  the  system  requests 
the  fifth  set  of  input  parameters  by  printing  out  the  phrase 
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INPUT  XO,  YO,  Rl,  R2,  SO:  (ELLIPSE  PARAMETERS) 


The  analyst  then  inputs  the  five  quantities 

XO  = longitude  of  center  of  ellipse  (in  degrees) 

YO  = latitude  of  center  of  ellipse  (in  degrees) 

Rl  = semi-major  axis  of  ellipse  (in  nautical  miles) 

R2  = semi -minor  axis  of  ellipse  (in  nautical  miles) 

SO  = orientation  of  ellipse  relative  to  due  north  (in  degrees). 

If  the  analyst  inputs  the  value  1=2  then  the  system  requests 
the  fifth  set  of  input  parameters  by  printing  out  the  phrase 
INPUT  N : (NUMBER  OF  POINTS  IN  BOUNDARY) 

The  analyst  then  inputs  the  number  of  points  N (N  < 30)  used  to  define  the 
boundary  of  his  selected  convex  confidence  region.  The  system  then 
requests  the  next  set  of  input  parameters  by  printing  out  the  phrase 
INPUT  X,  Y:  (LONGITUDE  AND  LATITUDE  OF  BOUNDARY  POINT) 

The  analyst  then  inputs  the  longitude  and  latitude  (in  degrees)  of  a 
boundary  point  of  his  selected  set.  This  last  operation  will  be 
automatically  repeated  N times. 

If  the  analyst  inputs  the  value  1=3,  then  the  system  requests 
the  fifth  set  of  input  parameters  by  printing  out  the  phrase : 

INPUT  X3,  Y3,  Bl,  B2:  (WEDGE  PARAMETERS) 

ft 

The  analyst  then  inputs  the  four  quantities 

XB  = longitude  of  vertex  of  wedge  (in  degrees) 

YB  = latitude  of  vertex  of  wedge  (in  degrees) 

Bl  = line  of  bearing  of  first  great  circle  of  wedge  (in  degrees) 

32  = line  of  bearing  of  second  great  circle  of  wedge  (in  degrees). 

The  lines  of  bearing  are  given  relative  to  due  north  with  Bl  < 32, 
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unless  the  wedge  spans  due  north. 

After  the  parameters  which  define  the  initial  confidence  region 
have  been  entered,  the  system  plots  on  the  graphics  terminal  a 30  point 
approximation  to  P and  then  computes  and  displays  a 30  point 
approximation  to  the  boundary  of  updated  confidence  region  . 

GRAPHIC  OUTPUT 

Two  examples  of  UPDAT  graphic  output  are  described  in  this  section. 

The  basic  program  update  operations  are  performed  in  spherical  earth 
coordinates.  The  coordinate  transformation  which  was  used  to  map  the 
spherical  earth  coordinates  into  the  plane  system  of  the  graphic  display 
is  called  a true  view  transformation  and  is  one  of  the  display  options 
provided  by  the  Graphic  Analysis  and  Correlation  Terminal  (GACT)  System. 

Figure  3 and  Figure  4 are  examples  of  updates  of  an  elliptical 

confidence  region.  In  both  cases  the  initial  confidence  region  of  the 

target  is  an  ellipse  with  center  (X0,  Y0)  = (0,  0),  semi-major  axis  of 

50  nautical  miles,  and  semi-minor  axis  of  25  nautical  miles,  and  which  is 

oriented  at  45°  relative  to  due  north.  In  the  first  example  (Figure  3) 

o o 

the  target  is  assumed  to  have  a heading  which  lies  between  0 and  30 
and  a speed  which  lies  between  20  and  30  knots. 

In  the  second  example,  (Figure  4)  the  target  is  assumed  to  have  a 
speed  which  lies  between  20  and  30  knots,  and  a heading  which  is  completely 
arbitrary.  In  both  examples  the  update  time  is  8 hours.  For  both  cases 
the  constraint  set  in  velocity  space  was  partitioned  into  a disjoint 
union  of  almost  convex  sets  and  a separate  update  operation  was  performed 
using  these  sets. 

The  total  updated  confidence  region  is  presented  as  a superposition 
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LAT,  LON:  0 0 

INPUT  XO,  YO,  Rl,  R2,  SO:  0 0 50  25  4 5 
INPUT  Al,  A2,  VI,  V2,  T .-0  30  20  30  8 


Fig.  3 — Example  of  updated  confidence 
region,  estricted  heading  bounds 
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LAT,  LON:  0 0 

INPUT  XO,  YO,  Rl,  R2,  SO  : 0 0 50  25  45 
INPUT  Al,  A2,  VI,  V2,  T : 0 360  20  30  8 


Fig.  4 — Example  of  updated  confidence 
region,  arbitrary  heading 
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1 


of  adjacent  updated  confidence  regions.  It  can  be  seen  that  the 
convexity  criterion  (see  Appendix  C)  could  have  been  relaxed  without 
losing  much  definition  in  the  updated  confidence  regions. 
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APPENDIX  A 


PROGRAM  LISTING 

This  appendix  contains  a listing  of  the  FORTRAN  version  of  the 
UPDAT  program.  The  subroutine  CONDRAW  used  in  the  graphic  version 
to  generate  a plot  of  the  polygonal  set  associated  with  an  array  of 
points  is  not  shown.  The  subroutine  which  computes  the  Time  View 
transformation  from  spherical  earth  coordinates  to  the  plane  graphic 
system  is  also  not  shown. 
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" GGRAM 


j 1 0 0 
0110 
0120 

00130 

' 0 14  0 

'10150 

0160 

00170 

00130 

0 0 1 9 0 
0 0200 

1 02 1 0 
.0220 
0 0230 
0 024  0 
0 0250 
00260 
0 0270 
■0  0230 
0 0290 
00300 
0 0 31  0 
00320 
: 0330 
0 0340 
i.i  035  0 
0 0360 
0 0370 
00380 
0 0390 
0 0400 
. 0410 
0 0420 
0 0430 
0 0440 
:'!  045  0 
■j  0460 
0 0465 
00470 
00430 
00490 
00500 
0 0510 
00520 
0 0530 
00540 


CONUP 


PROGRAM  coni  < INPUT.  OUTPUT.  TAPED 
COMMDNyBL  OK  1 y XE  <100).  YE  < 1 0 0> 

COMMON.-' BL  0K2/U  < 1 0 O')  . V < 1 0 0) 

C0MM0NyBL0K3y.XU  <1  00>  . YU  <100) 
COMMaNyBLOK4yXCi.YG.Rl  .P2.S0.V1  .V2.T.P 
CaMMONyBLOKSyXl  <'3000)  >Y1  <3000)  »Z<3000)  .TFK50) 
READ  <1 . ) ND.X0.Y0.P1 .P2.S0 
READ  <1 > > ND.A1.A2.V1.V2.T.P 
PRINT  10 

10  FORMAT <5X»26HB0UNDAPY  OF  UPDATED  PESIONyy) 

CALL  TEST  (At . A2> 

STOP 

END 

SUBROUTINE  TEST<A1.A2) 

CGMMON/BLOK4/XO.YO.P1 .R2.SO.V1 . V2.T.P 
DV 1 =V2— V 1 
B V 2 = V 2 ♦-  < A2  - A 1 > 

VUtAMAXI < DV1 » DV2) 

VL=AMIN1 < DV 1 » DV2) 

IF<V1.LE.  0.  01)  GO  TO  10 
IF (RlyT.LE. 0. 05*VL>  GO  TO  20 
IF  <VU.LE.  0.  05*P2-'T>  GO  TO  30 
CALL  PAPTCAl > A2.V1) 

RETURN 

10  CALL  IJPDA1  <A1  . A2) 

RETURN 

2 0 CALL  UPDA3  <A1 . A2> 

RETURN 

30  CALL  UPDA4CA1.A2) 

RETURN 

END 

SUBPOUT  I NE  PART  <A 1 . A2 . V 1 > 

E=0. 04 

AL=2^hC0S  < 1 . -EyVl) 

NT  = I NT  ( < A2-A 1 y AL>  + 1 
DAL=  (A2-A1 > yNT 
I F <NT , GT . 1 > GO  TO  10 
CALL  IJPDA1  <A1 . A2.) 

RETURN 

10  DO  50  1=1. NT 

A 1 =A 1 +DAL ♦ < I - 1 ) 

A2=A1 +BAL 

50  CALL  UPDA1CA1.A2) 

RETURN 

END 

SUBROUTINE  UPDA1 CA1.A2) 
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01020 

V CIO  =V  CK 1 ) +BV 1 ♦S I M CR1  > 

0 103  0 

GO  TO  10 

0 1 040 

13 

TFCKK2-IO  15*14.14 

01050 

14 

Jl'Jl+1 

01060 

U CIO  =V2*€0S  Cfll+Jl  ♦D31> 

01  070 

V CK>  =V2*SINCR1-kJ1  ♦DRl) 

01030 

GO  TD  10 

01090 

15 

IFCKK3-K)  17*16.16 

0 1100 

16 

U CK>  =U CK 1 > -DV3*C0S  CR2) 

OHIO 

V CK)  =V (Kl) -DV3»SIN<R2> 

01120 

GD  TO  10 

0 1130 

17 

I F CK'k'4— K ) 19.13.13 

01140 

13 

J2=J2+1 

0 1 1 5 0 

U CIO  =V1  ♦COS  CR2-.J2*DR2> 

0116  0 

V CIO  =V  1 ♦S I N CR2-.J2^DR2> 

01170 

GO  TO  :o 

01130 

19 

M=K  — 1 

01190 

RETURN 

01200 

END 

01210 

SUBROUTINE  BOPCT»M.N.L> 

01220 

COMMON 'BLOK  1 'XE < 1 00) .YE<100> 

01230 

CQMM0N'BL0K2'U Cl 00) .V  Cl 00) 

01240 

C0MMGN'BL0K3'XU  C 1 00> . YU  Cl  00) 

01260 

COMMON'BLOKS'Xl  C3000)  *Y1  C3000)  »ZC3000> 

> T R C 5 0 ) 

01265 

MN=M^N 

01270 

DD  10  1=1. N 

01230 

DO  10  .J= 1 * M 

01290 

K=M^CI-l>  +J 

01300 

X 1 CIO  =XE  c I > +T ♦( J CJ> 

01310 

1 0 

Y 1 CK)  = YE  C I)  +T*V  CJ> 

01320 

DO  15  1=1*30 

01330 

15 

Tfi  C I ) =TRN  C 0 . 2 0944*  C I — 1 :• ) 

01340 

DO  30  1=1*30 

0 1350 

I T=  1 

01360 

TH=0.20944*CI-1> 

01370 

IFCCTH.GT. 1.5703) . RND. CTH. LT. 4. 71 24) ) 

IT=-1 

01330 

DO  20  J=l*  MN 

01390 

20 

Z C .J)  = IT»CY  1 C.J -TR  C I > ♦;<  1 < J>  > 

0 1400 

25 

10=1 

01410 

M 1 =MN— 1 

01420 

DO  60  J= 1 . M 1 

01430 

J1=J+1 

01440 

I F CZ  C.J  1 > -Z  C I 0)  ) 5 0 > 5 0*60 

01450 

50 

I0-J1 

01460 

60 

CONTINUE 

01470 

XU  C I > =X  1 C I 0) 

01430 

30 

YU  C I ) =Y  1 C I 0) 

01490 

L =3  0 
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i 4 '5*  .5 
1494 

1 4'94 
1 498 

0 1 5 0 0 
0 1510 

02  0 0 0 
0201  0 
02  02  0 
02030 
02  040 
02050 
02060 
02070 
021  0 0 
02 1 1 0 
02120 
02130 
02 14  0 
02150 
02160 
02 17  0 
02180 
02 1 9 0 
022  0 0 
022 1 0 
02220 
02230 
02240 
02250 
02270 
08280 
02290 
02300 
02310 
02320 


DO  31  I = 1 * M 
PRINT  30.  1.1  >:'  I > . V < I > 

3 0 FORMRT  C2X ,F1 0. 3 » 3X » F 1 0.  3> 

31  CONTINUE 

RETURN 

END 

SUBROUTINE  UPHft3 <B1 . R2> 

COMMON/ BL0K2/U  <1  00>  » V <1  00> 

COMMON/  BL0K3/XU  < 1 00  ■*  r VIJ  < 1 0 0> 

COMMON/ BLOK 4/X 0» Y0.Ri.R2. SO . V 1 . V2 t T > P 
CRLL  BOS  Cfil » R2  > VI .V2.T.P.M) 

DO  10  I = 1 . M 
XU<I>=X0+U<I)  *T 
10  YU < I > =Y  0+V < I > ♦T 
DO  36  1=1 . M 
PR  I NT  35 . XU  t l > » YU < I > 

35  FORMRT  C5X » F 1 0 . 3 j 5X  > F 1 0 . 3) 

36  CONTINUE 

RETURN 

END 

SUBROUTINE  UPDR4(R1.R2> 

COMMON/BLOK 1 /XE < 1 0 U> . YE  <1 0 0> 

COMMON/ BL0K3/XU < 1 0 0> » YU  < 1 0 0> 
COMMON/BLOK4/X  0 , Y 0 . R 1 > P2 » S 0 . V 1 . V2 rT  tP 
L HL  L Bi  0 E ( X 0 * Y 0 » P 1 . P 2 > P * S 0 . N ) 

VX=0 . 5*  i'Vl  +V2)  ♦COS  ( 0 . 5 ♦ ■' R 1 +R2>  > 

VY=  0 . 5 ♦ < V 1 +V2> ♦S I N < 0 . 5 ♦ <B 1 +R2>  > 

DO  10  I = 1 . N 

xu«rn=xEa>+vx#T 

1 0 YU < I> =YE CI> +VY>T 
DO  36  I = 1 » N 
PRINT  35 » XU  <I'>  . YU  •:  I> 

35  FORMRT  CSX . F 1 0 . 3 . 5X . F 1 0 . 3> 

36  CONTINUE 

RETURN 

END 
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APPENDIX  3 


MATHEMATICAL  BASIS 

This  appendix  shows  that  the  update  algorithms  yields  the  updated 
confidence  region  with  an  accuracy  dependent  upon  the  number  of  points 
used  to  describe  the  boundaries  of  the  constraint  sets  p and  (/ 
in  position  and  velocity  space.  For  any  set  S?  let  3S  denote 
the  boundary  of  the  set  and  let  9(x,v)  = x + tv  denote  the  update 
map.  Then 

?u  = je(x,v;  |xeP,vel/]  = 9(p,i/.  . 

The  following  elementary  results  are  needed: 

Proposition  1:  If  the  sets  P and  1/  are  convex,  then  so  is  p^. 

Proposition  2:  The  boundary  of  P is  a subset  of  the 

u 

update  of  the  boundaries  of  P and  (/  . 

That  is,  3 P c 9(3  p,  3 \J) . 
u ' 

To  prove  Proposition  1,  let  ^ , ~J_2  e p • 

Then 

*1  = 5l  + %. 

= *2  + t-2 


•with 
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But 


\2Cj_  + (1~\)X2  S P 


\v  + (l-\)v2  e V. 


Thus 


>*1  * (1'^2  « Pu  ' 


To  prove  Proposition  2,  suppose  ^ = Xq  + tv^  s 3 Pu  . It  will  be 

shown  that  if  either  x / Bp  or  y i 5 1/ , a contradiction  is  reached. 

Suppose  ^ 3 p.  Then  there  exists  an  open  neighborhood  N..  about 

— o 

x With  N c p.  The  restriction  6 lN„  xv  of  9 to  the  set 
J x — o 

— o — o 

N x v is  a homeomorphism  and  so  mans  interior  noints  of  N X v into 

x -o 

interior  points.  This  implies  that  x^  + tv^  is  an  interior  point  of 

p which  is  a contradiction.  Likewise  the  assumption  that  v i a ti 

u * — o r ^ 
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leads  to  a contradiction.  Thus  x e 3 P and  v s . 

-o  -o 

The  computer  algorithms  which  executes  the  update  operation 

proceeds  as  follows : 

(a)  Approximate  the  boundaries  3P  and  3 \]  of  P and  V 

by  discrete  sets  Op  )Q  and  0 \J  ) consisting  of  30 

equally  spaced  points. 

(b)  Determine  the  set  9((3p)  , (3  1/)  ). 

o o 

(c)  Given  a finite  set  S = 9 ( (3  P ) , (3  V ) ) of  points, 
identify  those  points  which  are  boundary  points  of  the 
convex  hull  of  S. 

Propositions  1 and  2 assert  that  all  boundary  points  of  the  convex 

set  p are  contained  in  the  set  0(3  P,d  l/).  Since  0 ( (3  P)  , (31/  ) ) 
u 00 

approximates  the  set  9(3  p,3  1/),  step  c)  applied  to  9((c  p)  ,(3  l/)Q) 

yields  an  approximation  to  the  boundary  of  p . Step  (c)  is  discussed 

u 


in  Appendix 


which  follows. 


APPENDIX  C 


ALGORITHMS  USED 

This  appendix  contains  a general  description  of  two  of  the 
algorithms  used  in  the  update  transformation.  The  first  algorithm 
yields  the  boundary  points  of  the  convex  hull  associated  with  a finite 
set  of  points  and  may  have  some  general  interest.  The  second  algorithm 
which  is  discussed  decomposes  a non-convex  set  in  velocity  space  into 
a union  of  almost  convex  sets  and  is  more  specialised  in  its  application 
For  any  set  S in  the  plane  let  C , S)  denote  the  convex  hull  of 
5.  The  algorithm  used  to  identify  the  boundary  points  dC(S)  of  C ' S ) 
can  be  understood  by  referring  to  Figure  Cl. 


The  special  form  of  the  algorithm  which  has  been  programmed  yields 

pTT 

at  most  30  distinct  points  of  dc( S).  Let  0 = f^(k-l)  and  let  L. 

K 3'“'  A 

denote  a line  having  slope  m,^  = tan  e^.  Suppose  > 0.  The  set  of 

all  y-axis  intercepts  associated  with  the  set  of  lines  L (k  fixed) 

which  pass  through  all  the  point  of  S is  determined.  That  point 

(xk,yk)  of  S for  which  the  y-axis  intercept  z (Figure  C-l)  is 
"til 

minimized  is  the  k point  of  3c(s).  Now  increase  k by  one  and  repeat, 
til 

If  < 0,  the  k point  is  obtained  by  maximizing  the  set  of  all  y-axis 

intercepts.  Note  that  for  all  k,  9 never  takes  on  the  value  5 

K 2 

3TT 

or  Ihe  set  i(x1,y1),  ^X30,y30^  contains  at  most 

30  distinct  points  of  dc(S). 

An  almost  convex  set  in  velocity  space  can  be  defined  relative  to 
the  maximum  distance  <^max  between  that  set  and  its  convex  hull.  A 
set  can  be  said  to  be  almost  convex  up  to  a specified  level  E if 
^max  ~ FiSure  02  depicts  a constraint  set, 

V = jv|v1  < |v|  < v2,Al  < arg  (v)  < Ag  } 

in  velocity  space  which  is  very  non-convex. 
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Decomposition  of  a Non  Convex  Set  into  a Union  of  Almost  Convex  Sets 


a which  determines  the  size  of  the  sets  l/.  is  related  to 


